Visualization of Branch Points in PT-Symmetric Waveguides 
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The visualization of an exceptional point in a VT symmetric directional coupler(DC) is demon- 
strated. In such a system the exceptional point can be probed by varying only a single parameter. 
Using the Rayleigh-Schrodinger perturbation theory we prove that the spectrum of a VT symmetric 
Hamiltonian is real as long as the radius of convergence has not been reached. We also show how one 
can use a VT symmetric DC to measure the radius of convergence for non VT symmetric structures. 
For such systems the physical meaning of the rather mathematical term: radius of convergence, is 
exemplified. 



In the past several years, following the seminal pa- 
per by Bender and Boettcher non-hermitian VT - 
symmetric Hamiltonians have caught a lot of attention 
(see Q and references therein). Under certain conditions 
VT - symmetric Hamiltonians can have a completely real 
spectrum and thus can serve, under the appropriate in- 
ner products, as the Hamiltonians for unitary quantum 
systems Q. 

Recently, the realization of VT - symmetric " Hamilto- 
nians" has been studied using optical waveguides with 
complex refractive indices [J, [5j . The equivalence of the 
Maxwell and Schrodingcr equations in certain regimes 
provides a physical system in which the properties of 
VT -symmetric operators can be studied and exempli- 
fied. 

An extremely interesting property of VT - symmetric 
operators stems from the anti-linearity of the time sym- 
metry operator. Consider a VT - symmetric operator H, 
i.e., [VT, H] = 0. Due to the non-linearity of T one can- 
not in general choose simultaneous eigenfunctions of the 
operators VT and H. However, if an eigenvalue of the 
H is real then it's corresponding eigenfunction is a also 
an eigenfunction of the VT operator. This property has 
come to be known as exact / spontaneously-broken VT - 
symmetry. Exact VT - symmetry refers to the case when 
every eigenfunction of the VT symmetric operator is also 
an eigenfunction of the VT operator. In any other case 
the VT - symmetry is said to be broken. 

Usually, the transition between exact and 
spontaneously-broken VT symmetry can be controlled 
by a parameter in the Hamiltonian. This parameter 
serves as a measure of the non-hermiticity. An important 
class of VT - symmetric Hamiltonians are of the form: 
H(X) = H o + iXV. Where iJ (and V) are real and 
symmetric(anti-symmetric) with respect to parity so 
that [VT,H] = 0. When A = the Hamiltonian is 
hermitian and the entire spectrum is real. The spectrum 
remains real even when A ^ as long as A < A c . At 
this critical value and beyond, pairs of eigenvalues 
collide and become complex, see for example [6|. Bender 



et al. Q showed that the reality of the spectrum is 
explained by the real secular equations one can write 
for VT - symmetric matrices. These secular equations 
will depend on the non-hermiticity parameter and, 
consequently, yield either real or complex solutions. 
Delabaere et al. jgj showed for the one-parameter family 
of complex cubic oscillators that pairs of eigenvalues 
cross each other at Bender and Wu branch points. Dorcy 
et al. 0| after proving the reality of the spectrum for 
a family of VT -symmetric Hamiltonians, showed (lpj 
that at the point where the energy levels cross, i.e., the 
critical value of A, a super-symmetry is broken and not 
only the eigenvalues but also the eigenfunctions become 
the same. 

Consider the family of Hamiltonians given by: H = 
Hq + XV = Hq + i|A|V. The existence of a branch 
point in the Hamiltonian's spectrum determines the ra- 
dius of convergence of a series expansion of the energy 
in A. Friedland and one of us [ll| proved that for two 
real symmetric matrices Ho and V that do not com- 
mute there exists at least one branch point Xb p for which 
^rU=Ai, p = oo. Therefore, the expansion of the energy 
in powers of A converges only as long as |A| < |Ab p |. The 
most common situation is when the branch point is as- 
sociated with the coalescence of two eigenfunctions and 
the two corresponding eigenvalues. Such a point in the 
spectrum is often referred to as an exceptional point [l3| . 
Exceptional points in physical systems have been stud- 
ied, e.g., in 3, EH, 0, E3, E3]- Recently, exceptional 
points have been observed experimentally in microwave 
cavities In general, as stated in the theorem above, 
the value of the parameter A at which the branch point 
occurs is a complex number. This demands the control 
of the Hamiltonian by at least two parameters [l5| . From 
the evidence of branch points in previous studies of VT - 
symmetric Hamiltonians, e.g., [6j], it is plausible to as- 
sume that the branch point is located on the imaginary 
axis, i.e., Xb p = i\Xb p \. Hence, the two parameter depen- 
dence reduces to a dependence on a single real parameter. 
Therefore, the study of exceptional points in VT - sym- 
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metric systems is strongly simplified. 

We restrict the discussion to the cases for which | Xb p | ^ 
0. Thus for A = i|A| and Xb p = i\\b P \ in the vicinity of 
an exceptional point one can write, 

E(X) = E»P±D^(X-\ bp )(\-\* bp ) (1) 

= ^±^|A bp | 2 -|A| 2 . (2) 

If |A| > \\b P \ one gets a pair of complex conjugate eigen- 
values E and E* , whereas in the case for which |A| < \Xb P \ 
the system has two real eigenvalues. Therefore, Xb p cor- 
responds to the critical value, A c , at which the transition 
between exact and spontaneously broken VT symmetry 
occurs. An interesting consequence of the above analysis 
is that if one treats the potential AV^ as a perturbation 
the Rayleigh-Schrodinger perturbation expansion would 
converge as long as |A| < |A& P |. One would then expect 
that the reality of the spectrum should be apparent from 
such an expansion. We shall prove this in the following. 

Consider the following time-independent Schrodinger 
equation: 

(H + XV)* j (X)=E j (X)* j (\), (3) 

Where again, iJo(and V) are real and symmetric (anti- 
symmetric) with respect to parity and are both hermitian 
operators. Assuming that |A| < |A& P |, one can expand 
the eigenvalues and cigcnfunctions in a convergent power 
series in powers of A: 

oo 

E ] {\X\<\K\) = E^^' ( 4 ) 

*j(|A| < I VI) = f> n *J n \ (5) 

where -Ey ) ; n = 0, 1, 2.. are the real energy (wave- 

function) correction terms of the Rayleigh-Schrodinger 
perturbation expansion. The (2n+l)- rule stated by 
Wigner [l2| implies that, 

E 3 = ( X ^\H Q + XV\ X ^)+0(X 2n+2 ) (6) 

where, 
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x5-%) = E A ^f (?) 

fe=0 

and (x^ IXj™"*) = 1- Therefore, following the (2n+l)- 
rule, 

E f n+1) = (^ n) \V\^ n) ) (8) 

where the n-th order correction to the exact eigenfunction 
$ j is the solution of the following equation, 

^;\x) = G' (Ef)V(x)\^) (9) 

ri-1 

fc=i 



andG^f ) =J2 q ^{x\4 0) )(4 0) \(Ef -E^)-K The 
parity symmetry of the zeroth-order Hamiltonian ensures 
that its eigenfunction ipq ^ (x) has either even or odd par- 
ity. From the hermiticity of Ho the eigenvalues of the 
zeroth-order Hamiltonian are real. Therefore, the n-th 
order energy correction term, {E^}, is also real. As 
can be verified by substitution, the definite parity of the 
ipq ^ (x) together with the antisymmetric parity of the 
perturbing potential, i.e., V(x) = — V(— x), generate a 
definite parity n-th correction to the wavefunction. That 
is, 

^ n \x) = (-iy +n - 1 ^ n \-x), (io) 

for: j = 1, 2, . . . and n = 0, 1, . . .. One can now immedi- 
ately conclude that, 

£f n+1) =0. (11) 

Therefore, 

oo 

^■(|A|<|A, p |) = E A2 "^r ) - ( 12 ) 
Consequently for A — i\X\ the series 

oo 

e j(\M < I'M) = ^{-l) n \X\ 2n Ef n) (13) 

converges to real values. 

The above analysis shows that the transition between 
exact and spontaneously broken VT - symmetry occurs 
at a branch point where the non-hermiticity parameter 
(|A| in our case) reaches the radius of convergence of the 
Rayleigh-Schrodinger perturbation expansion. As we will 
show below this phenomenon can be observed experimen- 
tally. 

The measurement of the radius of convergence of a 
Rayleigh-Schrodinger series expansion will be demon- 
strated in a VT - symmetric waveguide configuration. A 
VT - symmetric waveguide can be easily realized with 
a symmetric index guiding profile and an antisymmetric 
gain/loss profile, i.e., n(x) = n*(—x) 0. We consider two 
coupled planar waveguides depicted in Fig. [1] for which 
the refractive index varies only in the x direction. The 
direction of propagation in the waveguides is taken to be 
the z axis. The wave equation for the transverse-electric 
(TE) modes then reads: 

[^ + k 2 n{xf^£ v {x)^p 2 £ y {x) 1 (14) 

where the y component of the electric field is given 
by: E v {x,z,t) = £ y {x)e l ^ fiz \ k = 2tt/A, and A is 
the vacuum wavelength. Clearly, the wave equation for 
the y component of the electric field, i.e., Eq. (TT4")) . is 
analogous to the one-dimensional Schrodinger equation: 
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FIG. 1: (color online) A VT -symmetric directional coupler. 
The structure consists of two coupled slab waveguides. The 
real and imaginary part of the refractive index is also por- 
trayed. The refractive index only varies in the x direction. 
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tial as V{x) — — \k 2 r? and the energy as: E = — i/3 2 
As shown in Fig. [T] we couple between one gain-guiding 
waveguide (positive imaginary part of the refractive in- 
dex) and one loss-guiding waveguide (negative imaginary 
part of the refractive index) [20| in order to create the 
VT - symmetric structure. For simplicity we take the 
separation between the two coupled waveguides to be 
the same as the waveguides' width, i.e., 2a. The effect 
of changing the separation of two coupled VT - symmet- 
ric waveguides was studied in [5(, and showed that as 
in regular (non VT - symmetric) directional couplers the 
coupling length is an increasing monotonic function of 
the separation. The beat time period usually used to de- 
scribe quantum beating systems , e.g., [2l| is exchanged 



in the study of waveguides by a beat length: L = 2ir/ A/3, 
where A/3 is the difference between the propagation con- 
stants of the two modes. As described below a different 
method of controlling the beat length is by changing the 
non-hermiticity of the potential which can be controlled 
by the gain(loss) coefficient: Aa. 

In order to illustrate the control of the beat length us- 
ing the non-hermiticity parameter we choose the follow- 
ing parameters for the waveguide structure shown in Fig. 
[TJ The background index is taken to be: no = 3.3, the 
vacuum wavelength: A = 1.55/im, the real index differ- 
ence between the waveguides and the background mate- 
rial: An = 10~ 3 , and the separation between the waveg- 
uides which equals the waveguides' width: 2a = 5/jm. 
The parameters are chosen such that each waveguide con- 
tains only a single guided mode before we couple them. 
The coupled guided modes are calculated by diagonal- 
izing the matrix representation of Eq. [14] in a sine ba- 
sis. The " Hamiltonian" matrix is non-hermitian and one 
needs to take care when normalizing the eigenvectors. We 
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FIG. 2: (color online) The two trapped modes of the wave- 
guide depicted in Fig. [T]as a function of the non-hermiticity 
parameter strength. The eigenmodes approach each other on 
the real axis as Aa increases until a critical value of Aa c ~ 8.4 
is reached. At the critical value one finds a branch (ex- 
ceptional) point where the two modes coalesce. Beyond the 
branch point the directional coupler sustains one gain guiding 
mode and one loss guiding mode. 



choose to normalize our eigenvectors according to the so- 
called c-product |22|, i.e., (£„|£ TO ) = (££|£m) = S n ,m- 

The coupled waveguides support two guided modes. 
The propagation constants of the two modes are plotted 
in Fig. [2] as a function of the non-hermiticity parame- 
ter. Increasing Aa causes the propagation constants of 
the two modes to move towards each other up to a criti- 
cal point, i.e., the branch point. Beyond the branch point 
the propagation constants become complex conjugates of 
one another. As long as VT - symmetry remains exact, 
i.e., Aa < Aa c ~ 8.4, the two guided modes can be 
classified according to the parity of the real part of the 
transverse electric field just as in the non VT - symmetric 
case. The critical value of Aa corresponds to a branch 
point, i.e., an exceptional point, where the two modes 
coalesce. At the branch point both the propagation con- 
stants and the corresponding electric field become equal. 
One can, therefore, study the exceptional point in a VT - 
symmetric waveguide by varying only a single parameter. 
Past the critical value of Aa, the waveguides support 
one gain-guiding mode and one loss-guiding mode. The 
transverse field past the branch point no longer retains 
the symmetry properties of the VT operator, but rather 
each of the two modes becomes localized in one of the 
waveguides. 

The progression of the propagation constants on the 
real axis towards the branch point can be visualized ex- 
perimentally by observing the beat length (time, in the 
corresponding quantum mechanical problem) of the sum 
field for the VT symmetric waveguide. Fig. [3] displays 
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FIG. 3: (color online) The power distribution for a propagat- 
ing sum field consisting of the two guided modes, see Eq. [15] 
The propagating light is shown for three values of Aa. As can 
be readily observed, the beat length (analogous to the beat 
time period in quantum mechanics) increases as the value of 
Aa approaches the critical value. See the text for further 
details. 



the power distribution, 



\E y ( X ,Z)\ 2 = 



1 
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(15) 



for three values of Aa. As the value of Aa approaches 
the critical value the beat length increases. This is a di- 
rect observation of the movement of the propagation con- 
stants towards each other on the real axis. Near the crit- 
ical value of Aa, i.e., the exceptional point, the sum field 
no longer oscillates between the waveguides but rather 
travels in both waveguides simultaneously. Therefore, by 
tuning the value of Aa one can visualize the movement 
of the eigenmodes towards the branch point or away from 
it, as seen in Fig. [3] 

The effect of the non-hermiticity on the beat length is 
clearly reflected in Fig. [3J However, this visualization of 
the branch point is not the only information one can get 
from such an experiment. By finding what is the criti- 
cal value of Aa which corresponds to the branch point 
one can infer the radius of convergence for any symmet- 
ric waveguide structure with an antisymmetric pertur- 
bation. Therefore, while the visualization is done on a 
VT symmetric waveguide, the critical value of Aa corre- 
sponds to the maximum value of an added antisymmetric 
index profile which can still be treated within perturba- 
tion theory. 

In conclusion, a proof of the reality of the spectrum of 
VT - symmetric Hamiltonians has been given using the 
Rayleigh-Schrodinger perturbation theory. We showed 
that the transition between exact and spontaneously bro- 



ken VT - symmetry occurs at a branch point and the 
transition can be visualized in a VT symmetric waveg- 
uide structure. Such a system has many advantages for 
the study of exceptional points as the movement towards 
the branch point is controlled by a single parameter. We 
also show that by such a measurement one can find the 
radius of convergence for the family of waveguides where 
a symmetric index profile is supplemented with an anti- 
symmetric index profile. This antisymmetric profile need 
not be imaginary. The knowledge of the radius of con- 
vergence is important if one wishes to treat the added 
antisymmetric index as a perturbation. 
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